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I. INTRODUCTION 

Since the seminal work of Sutherland and Shastry jl],|| , 
the ground state properties in quantum systems with adi- 
abatic change of twisted boundary conditions have at- 
tracted much current interest. Although the period of 
the full spectrum should be 2n as a function of the twist 
angle, each eigenstate has a period larger than 2n in gen- 
eral, if we trace the flow of the initial state adiabatically. 
This period provides us various information for quantum 
spin systems, correlated electron systems, etc. For ex- 
ample, it can be a sensitive probe to examine whether 
particles form a bound state like Cooper pair d^]. 

One of the remarkable applications of this method is 
to observe (fractional) exclusion statistics || directly in 
the spectral flow of the ground state for one-dimensional 
quantum systems For this purpose, it is crucial to 
sweep away all irrelevant interactions from the model, 
because the period is quite sensitive to irrelevant inter- 
actions, leading to the period 47T if they exist and 
hiding the nature of exclusion statistics. Among oth- 
ers, it is known that quantum many-body systems with 
1 /r 2 interaction [|]]|] provide us with an ideal situation to 
observe exclusion statistics, i.e., these systems are com- 
pletely free from irrelevant interactions. Motivated by 
this, we have recently studied the spectral flow in the 
Haldane-Shastry spin model with 1/r 2 exchange by gen- 
eralizing the model to include twisted boundary condi- 
tions Q. We have indeed found that the period of the 
spectral flow is controlled by the statistical interaction in 
exclusion statistics. 

In order to investigate exclusion statistics for corre- 
lated electron systems, we wish to ask a basic question, 
i.e. what happens for the spectral flow when holes are 
doped into the Haldane-Shastry spin model. To address 



this question, we investigate in this paper the exact spec- 
tral flow in the supersymmetric t-J model with 1/r 2 hop- 
ping and exchange, which was exactly solved for the first 
time by Kuramoto and Yokoyama with periodic bound- 
ary conditions 10 -ill. We extend the model to include 
the twisted boundary conditions, and then discuss spec- 
tral properties of the model as a function of the twist 
angle. 

The paper is organized as follows. In the next sec- 
tion, we briefly mention how to generalize the solvable 
supersymmetric t- J model of 1/r 2 interaction to be com- 
patible with twisted boundary conditions. In section III, 
we then derive the exact eigenstates of the model by ex- 
ploiting a Jastrow-type ansatz. In section IV, it is shown 
that the spectrum thus obtained can be given alterna- 
tively by the asymptotic Bethe ansatz (ABA). It is then 
demonstrated that the spectral flow of the model is com- 
pletely described in terms of the motif in the ABA. We 
find that once holes are doped into the Haldane-Shastry 
spin chain, the nature of fractional statistics can not be 
observed in the period of the spectral flow, though ele- 
mentary excitations are still characterized by exclusion 
statistics. We claim that this is related to the breaking 
of the commensurability condition in the motif picture. 
Section V is devoted to summary. 



II. MODEL HAMILTONIAN 



We consider the supersymmetric t-J model with 1/r 2 
interaction The model Hamiltonian is given by 
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where T 5 is the projection operator to exclude the dou- 
ble occupancies at each site, i and j denote the indices 
for sites — 1, 2, • • • N), and a is a spin index labeled 
by a =1 and j. Note that the summation with respect 
to I is introduced to treat the system in a ring geome- 
try Kuramoto and Yokoyama jHj] have shown that 
this model can be solved exactly with periodic bound- 
ary conditions, in which the summation over I results in 
the sine-inverse-square interaction. This model, which 
includes the Haldane-Shastry spin chain at half filling, 
provides us with a paradigm of ideal exclusion statistics 

IS 

Preserving the integrability, we wish to treat the model 
with twisted boundary conditions. Note that this type 
of generalization is not straightforward for the present 
model because of the long-range nature of hopping and 
interaction. As was the case for the spin chain ||, it 
turns out that the summation over I in eq.(|l|) after im- 
posing twisted boundary conditions plays a crucial role 
to preserve the integrability. The boundary conditions 
we impose in this paper are 
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In what follows, the angle <j) a in unit of 2tt defined above 
is referred to as the twist angle. Note that we can impose 
twisted boundary conditions independently on the spin 
and the charge degrees of freedom with two kinds of twist 
angles. To diagonalize the Hamiltonian, it is convenient 
to treat the system after a gauge transformation, 
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as was done previously . Hence the transformed oper- 
ators satisfy the periodic boundary condition. 

In the remainder of the paper, we restrict ourselves to 
the case of rational twist angles, 
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where p a and q a are integers. This restriction enables 
us to carry out the summation over / in eq.(|l|), and the 
resultant Hamiltonian is well defined, 
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where 
is 



V, (5) 



and the effective coupling constant 
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for <j) = p/q. The effect of twisted boundary conditions 
is now incorporated in the effective coupling (||), and we 
can solve (|) with periodic boundary conditions. Putting 
p = 0, we have Jo{n) — sin -2 reproducing the peri- 
odic model [[To|| . 

In order to obtain the exact eigenstates of the Hamilto- 
nian (||), we take the fully polarized state \P) = \ 11 ■ ■ ■ 1) 
as a reference state, and introduce the fermionic hole op- 
erator hj by |l(J 

Ci| = h\, c] t = hi, (7) 
and the bosonic spin operator bj by 

Si = 4, S+=b tl S?=^--blb i , (8) 

where = 1 — h\hi. The Hamiltonian now reads 

H= (^)VT(0T)+Wl) + ^) + #int+e), (9) 
where 
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Here and in what follows, we omit V for simplicity, 
though we regard the Hamiltonian as always accompa- 
nying this operator. 



III. EXACT SOLUTION FOR TWISTED 
BOUNDARIES 

In this section, we obtain the exact eigenstates of d) 
by exploiting a Jastrow-type ansatz for the correspond- 
ing wave functions ||s]-|Ic| . Let us first divide the total set 
of sites {1, 2, , • • ■ , N} into three subsets {a; Q }, {yi} and 
{u a }, which denote, respectively, the position of down- 
spins (bosons), holes (fermions) and up-spins ||lO,|l2|. We 
use Greek and Latin indices for bosons and fermions, re- 
spectively, the numbers of which are denoted by M and 
Q. Then, we have S z = ^f - - M. 



2 



As was found in |l0| , electronic wave functions of Jas- 
trow type are the exact eigenfunctions for the t-J model 
with periodic boundary conditions. We can also expect 
such wave functions to be still eigenstates even for twisted 
boundaries (Q) if we treat the twisted systems after the 
gauge transformation (|3|), as was indeed demonstrated 
for the Haldane-Shastry spin model These observa- 
tions motivate us to consider a Jastrow-ansatz state as 
an eigenstate for the twisted t-J model, 
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where the wave function ip is defined by 

ip({x a }, { yi y, j s , j h ) = n n z j -y 



(11) 



a 



I 



x JJ d(x a - xfjf Y[ d{yi - y m ) \\ d(x a - yi ). (12) 
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Here, d(n) — sin(7m/AT), and we have introduced 

z = e 2 "/ N . (13) 

Spin and charge currents are defined by J s = Q/2 mod 
1 and J c — (Q + M — l)/2 mod 1. Note that the wave 
function (|l^) is the gauge-transformed one whose basis is 
constructed by the periodic operators b Xa+N = b\. a and 

hy [+N — b^ t (see (pi])). The original wave function, based 
on the twisted operators, is obtained via the transforma- 
tion (|), 

ipoT g ({x a },{yi}; J s , Jh) 

= n ^ n z * m ■ j « j ^ (1 4 ) 

a I 

We will show below that the wave function is indeed an 
eigenstate of the present Hamiltonian. For this purpose, 
we first calculate the actions of hopping terms on the 
wave function ( |l2"| ) in several steps, following similar cal- 
culations to Refs. jlO| and (l^] except that hopping terms 
are now given by a complicated form depending on cj) a 
in a nontrivial way. 



A. Action of T T ( 



The action of Tf on the wave function is 
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Here we have introduced the same notations as those in 

Pi. 



B. Action of T 



Note that exchanges pairs of {x a } and {y{\, and 
is difficult to treat directly [To|JT2]| . However, if we in- 
troduce a unitary operator U generating the 7r-rotation 

around ir-axis; U = ■ e m ^ S j +S i can be expressed 

by 
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Furthermore, from the identity proved by Wang-Liu- 
Coleman Jl|, 

JJ d 2 (x a - Xp) \\_d{x a - y{) = 
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where A is a constant independent of the coordinate, we 
find a generalized Wang-Liu-Coleman's theorem includ- 
ing the twist angle, 

Ti(<fiiM{x},{y};J s ,J h ) _ 
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Thus, the action of Tj. can be evaluated in a similar way 
to Tf as 
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where J| = — J s + 7V/2 and is defined as eq.(fi6| 
replacing y m by w Q in eq. ([!]]). 



C. Action of T s {4>s) 

The action of T s on the wave function is 
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Here we denote z a = z Xa 



D. Eigenvalues 



There are various many-body terms in eqs.(|l 



and (22). However, more than three-body terms turn out 
to vanish similarly to the case of the periodic model. Key 
identities to prove this are 
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where J is an integer < J < AT, s and i are non-negative 
integers in the region < s + 1 < min( J + <f>, N — J — <f>). 
The function e{k) for any rational k is defined by 



e(k) = (2[fc] - N + 1) k - [k] {[k] + 1) 



(26) 



where [a] denotes the Gauss symbol, denoting the max- 
imum integer which does not exceed a. It is quite char- 
acteristic in twisted 1/r 2 models that such a function 
including the Gauss symbol enters the expression. We 
will see later that e(k) is nothing but a single-particle 
energy as a function of the momentum k. 

Substituting the above results to eqs.([l5|); (|2ll ) and 
(^||), we have 
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First, Wq is the constant term 



Wo = -M(N 2 - 1) + -M(M 2 - 1) + -QM(2M - 1) 
3 3 2 

+ ^Q(N 2 - 1) + 1q(Q 2 - 1) + ±Q 2 (N - Q) 



+2Me(J s + <f> s ) + 2Qe(J h + T ) 
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where e( J + (f>) is defined by 
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= { l{ £ {J + 4>~\) (29) 

+e( J + 4> + 5) — 5), for half-integer J. 



Next, W2 is the two-body term, evaluated as 
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which should be canceled out by one of the interaction 
terms in eq. (|To|) . Finally, the three-body term W3 is 

W 3 = - Y cot cot + \ E cot ® la cot ® ma 

= ~2J ^(^a -yO 

a,; 

+~Q(M + 1) + iQ(7V 2 - 1) 
2 

-l Q 2( N -M) + ~Q(Q 2 -l). (31) 

The three-body term thus reduces to the two-body term, 
and is also canceled out by the other interaction terms 
in eq.(|l^). Collecting constant terms, we finally end up 
with the eigenvalue of the Hamiltonian 
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with 



e = ^Q(N 2 - 1), 
E = ^M(M 2 - 1) 

+ \Q{Q 2 - 1) + \Q{M + Q)(2M - Q) 
+ 2Me(J 8 + T - ^) 
+ 2Qe(J h + (f> r ) 

+ 2Q(e(J l + ( f >l ) + — 



(33) 



Here we have divided the energy into two parts e and 
E for later convenience, because E can be reproduced 
by the asymptotic Bethe ansatz (see eq.(|39|)), whereas e 
simply corresponds to the chemical-potential term. The 
currents are subject to the constraints 
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in order for the formula (|25|) to be applicable. Note that 
when <j)^ — (pi — 0, e(J) reduces to the well-known form 
s(J) = J (J — N), and the spectrum obtained here coin- 
cides with that found in JT^ ]. 

The expression ( |33"| ) gives the exact ground state en- 
ergy with twisted boundary conditions. However, in or- 
der to trace the flow of the initial ground state correctly, 
we need more detailed information to specify the states, 
which may be supplied by the motif picture described 
below |1. 



IV. SPECTRAL FLOW 



010- ■ -010, where the arrow means 5(f> — 1, or in other 
words the ring is threaded by a unit flux. Thus the period 
of the ground state turns out to be 2, which directly re- 
flects the statistical interaction g — 2 of exclusion statis- 
tics f§. 

Keeping the above results in mind, we now show 
how the ABA equation for twisted boundaries is gen- 
eralized when holes are doped into the Haldane-Shastry 
model. To this end, we exploit the idea of the nested 
Bethe ansatz with graded symmetry (supersymmetry) 
]i~6[ , which leads to the Bethe equation for the super- 
symmetric t-J model with twisted boundary conditions, 



M+Q 



r(2) 



T^sgnfe -TO;) + s S n ( fc M-M> (37) 

i=l u=l 
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- = o s S n (™' ~ M> (38) 



A. Asymptotic Bethe Ansatz 

Before discussing the spectral flow as a function of the 
twist angle, we first show that the above exact spec- 
trum can be reproduced by the asymptotic Bethe ansatz 
(ABA), which naturally leads us to introduce the motif 
picture. To obtain the ABA equation for the supersym- 
metric t-J model, we recall that the present model is 
reduced to the Haldane-Shastry model at Q — 0, and the 
spectral flow in this case HQ can be described by the 
ABA equation, 



kfj, — I ^ 
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We then have a simple relation — + <p s , in terms 
of the rapidities fc's defined at 4> s = 0. Namely, the sole 
effect due to 4> s is to shift k's uniformly. To observe how 
the spectral flow occurs, it is convenient to introduce the 
description by the motif fl5|] , which is briefly outlined 
here in the case of half filling. The essential point of 
the motif in this case is that one can describe the effect 
of two-body phase shifts in (^5|) by arranging and 1 
for a given configuration of rapidities. For example, the 
ground state of the Haldane-Shastry model is a singlet, 
denoted by the motif 



010- ■ -1010101 •••010, 



which means that the spacing of the occupied momenta, 
denoted by l's, is enlarged twice as large as that for free 
fermions. From this motif we can identify the statistical 
interaction g = 2 (g = 1 corresponds to the free fermion 
case). The spectral flow of this state can be described 
by the flow of the motif, 010- • -010 -> 001 • • ■ 101 -> 



where \i — 1, 2, • • • , M+Q and i = 1, 2, • • • , Q. The newly 
introduced rapidity rhi is concerned with the charge de- 
grees of freedom. The energy E in eq.(B2h is given by 



E 



M+Q 

l"=l 



(39) 



where the single-particle energy e{k) is defined by eq.(|2f 

Note that the twist angles, cf>-\ and (fii, have been 
introduced in the above equations so as to be consis- 
tent with a nested procedure in the Bethe ansatz [ p^[ . 
Both_of the rapidities k and m are defined in the range 
< k, m < N, which include the effect of the finite twist 
angle. The corresponding quantum numbers 1^ and 1^ 

should satisfy = M/2 mod 1 and /f 3 = (M + Q)/2 
mod 1. 

Here we should mention the precise meaning of eq. J3S| ) 
which has been formally deduced by the nested Bethe 
ansatz, because it does not seem to hold for a fractional 
value of (f>i at a first glance. This subtle problem comes 
from the phase shift with the step-wise sgn(fc) function. 
We briefly summarize how to treat this equation cor- 
rectly. Consider, for example, the following configuration 
at </> T = <j>i = 0, 



• • • < feu-i < rrii < kn < 



(40) 



(36) If we put 



1, it is seen from eq. 



that I 



(2) 



(2) 

changes into D +1 and the above configuration should 
be changed to 



< k„ 



-i < k^ < m l 



(41) 



Namely, m, exchanges the position with k^ sitting on its 
right. Therefore, fractional cj>i between and 1 should 
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interpolate these two configurations smoothly. This is 
naturally realized if we introduce an infinitesimal width 
r\ in the step-wise phase shift and then take the limit of 
r\ — > 0. Based on this observation, we can separate the 

l.h.S. Of eq.© as 4 2) + 0; EE (if > + [0 X ]) + (0 X - [0;]), 

and the fractional part 0| — [0J can be absorbed into 
sgn(rni — kft) EE sgn(0) in the case of TOj = k ^ for 
< 0| < 1. Namely, the fractional portion of 0j can 
be incorporated into sgn(0) by taking the appropriate 
limit mentioned above ]l9| . As a consequence, we have 
from eq. (|37|) 



1 



sgn(mi 



(42) 



for < <fii < 1. Note that the initial order m 
01 = indeed changes into k^ < to, at 0^ = 1. Substi- 
tuting the above equation into eq. (|3^) , we end up with 
the simple results, 



kfj, — 



' 0T; 



for ■ • • < to, < kfi < 



for • • • < k,. 



-1 < fc p < 



(43) 



for < (f>i < 1. This equation implies that fc with (with- 
out) to in its left neighbor describes the charge (spin) 
degrees of freedom One can easily confirm that the 
above ABA equations indeed reproduce the exact energy 
( |33| ) obtained in sect. Ill by choosing the suitable quan- 
tum numbers 1^ and 1^ . 



B. Spectral Flow in Metallic Phase 

Having noticed that the exact solution obtained in 
sec. Ill can be reproduced by the ABA, we are now ready 
to discuss the spectral Sow by comparing the numeri- 
cal diagonalization results with the motif picture in the 
ABA. The numerical results for small systems can be 
complementary to the exact results in sect. Ill, because 
the analytical method there can supply only a special se- 
ries of the eigenstates. In Fig.l, we have shown the exact 
spectral flow calculated numerically for two holes in the 
system with 6 sites (i.e. 4 electrons). In this system, the 
ground state is doubly degenerate concerning the spin 
degrees of freedom. For instance, two degenerate ground 
states in Fig.|l| at 0f = <f>i = are classified by the motifs 



Ml 

M2 



OlololOlO, 
OlOlololO, 



and the first excited state by 

M3 : OlolOlolO. 



(44) 



(45) 



the position of to, in eqs.(|37|) and (pq), namely, it spec- 
ifies a doped hole in the motif. When boundaries are 
twisted, the pattern behavior in the spectral flow explic- 
itly depends both on the configurations of A:'s and to's, 
because the motions of k's arc different from those of 
to's as a function of the twist angle. These behaviors 
are correctly described by the flow of 1 and o in the 
motif. Since we can change two independent twist an- 
gles in the doped case, we have examined the spectral 
flow for two typical cases in Figjl]. One is the case where 
0T = (fii = 0, which will be referred to as Type I. In this 
case, the twisted boundary condition is imposed only on 
the charge degree of freedom. The other is the case where 
01 = 0, 4>i = 0, which will be referred to as Type II. It is 
seen from eq. (^3|) that the twisted boundary condition in 
Type II is equally imposed both on the charge and spin 
degrees of freedom. The motif Ml, for example, behaves 



OlololOlO -> OlOlololO for Type I, 
OlololOlO -> OOlololOl for Type II, 



(46) 



when we increase the twist angle 
corresponding spectral flow reads 



= 1. Therefore the 



Ml : al 
M2 : a2 
M3 : bl - 



• a2 
■ a3 
b2 - 



> a3 

► al 
63- 



> a4 (= al), 
♦ a2, 
61, 



for Type I shown in Fig. 1(a), and 



Ml : cl -> c2 
M2 : c3 -» c4 



c3 
c5 



c4 
c6 



c5 
cl 



c6 
c2 



cl, 
c3, 



(47) 



(48) 



for Type II (M3 is the same) shown in Fig. 1(b). Note 
that — + means (50 = 1. In this way we can trace the 
natural spectral flow correctly following the flow of the 
motif. This in turn demonstrates that the spectral flow in 
the supersymmetric t- J model with 1/r 2 interaction com- 
pletely fits in with the motif picture in the ABA even in 
the metallic case. We note that the spectral flows Ml 
and M2 for the ground states (Fig.|l|) agree with the ex- 
act results obtained in sect. III. 

Based on these results for small systems, we now wish 
to generalize our statement to more generic cases, and 
deduce how the spectral flow behaves when any number 
of holes are doped into the Haldane-Shastry model with 
TV-sites. We first note that the doped holes energetically 
prefer to occupy the positions in the center of the motif. 
For example, if we dope one or two holes, the correspond- 
ing motif should be 



010---1010ol01---010, 
010---101olol01---010. 



(49) 



In order to express the motif in the doped case, we have 
introduced o in addition to and 1, where o denotes 



It is easily seen from these motifs that the period of the 
flow becomes macroscopic, i.e. N both for Types I and 
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II. Note that this is also the case for any doping rates 
away from half filling, namely, the period generally be- 
comes TV both for Type I and Type II. The period N for 
Type I (the charge sector) is naturally understood by the 
analogy to the free fermion case. On the other hand, the 
change of the period in Type II is remarkable, because it 
is given by 2 at half filling, reflecting fractional exclusion 
statistics with statistical interaction g = 2. Therefore 
once holes are doped, characteristic properties of ideal 
exclusion statistics for the spin sector are hidden in the 
spectral flow. This does not imply, however, that the na- 
ture of exclusion statistics is spoiled by doping, because 
various characteristic properties such as correlation func- 
tions can be still interpreted in terms of exclusion statis- 
tics A reason why we cannot observe exclusion 
statistics in the spectral flow is related to breaking of 
the "commensurability" condition. For example, exclu- 
sion statistics with g — 2 is commensurate at half filling 
in the sense that the corresponding motif for the ground 
state fits in with the underlying lattice (see eq.(36)). In 
contrast, the motif for the metallic phase, for example, 
( [l9| ) does not meet the commensurability condition, hid- 
ing the g = 2 statistics in the flow. Therefore, we can 
say that the commensurability condition is essential to 
observe the ideal exclusion statistics in the spectral flow. 

In this connection, it is instructive to think of a specific 
motif for an excited state in the case of N/3 hole-doping, 
which has a uniform configuration such that 



01ol01ol0---01ol0. 



(50) 



This motif satisfies the commensurability condition, and 
the resultant period of the state turns out to be 3. Al- 
though this kind of a special excitation may not be in- 
teresting for the present model, we wish to note that the 
same motif describes the "ground state" for the SU(3) 
Haldane-Shastry spin chain. Therefore, we can predict 
that the ground state of the latter model, or more gener- 
ally, SU(^) spin model can have a period v in the spectral 
flow. 

The above characteristic properties for the 1/r 2 model 
are to be compared with ordinary interacting electron 
systems, for which the ideal statistics may not be real- 
ized due to perturbations which deviate the system from 
the idealized situation. In such cases, it is difficult to ob- 
serve exclusion statistics in the spectral flow even if the 
commensurability conditions are met, because irrelevant 
perturbations may dominate such properties. 



V. SUMMARY 

We have investigated the spectral flow in the super- 
symmetric t- J model with 1/r 2 interaction by analyzing 
the exact spectrum with twisted boundary conditions. 
At half filling, exclusion statistics with statistical inter- 
action 5 = 2 can be observed in the period of the spectral 



flow. Away from half filling, two kinds of spectral flows 
for the charge and spin sectors appear. It has been found 
that the nature of the spectral flows nicely fits in with 
the motif picture even in the metallic phase. We have 
also addressed the question whether one can observe frac- 
tional exclusion statistics explicitly in the spectral flow 
in the metallic phase. By analyzing the exact spectrum 
in terms of the motif picture, we found that the period 
becomes N once any number of holes are doped, hiding 
characteristic properties of exclusion statistics with g = 2 
for the spin sector. This drastic change has been shown 
to be related to breaking of the commensurability condi- 
tion. Therefore, even in the idealized model for exclusion 
statistics, which is free from any irrelevant perturbations, 
it is crucial to satisfy the commensurability condition to 
observe the statistical interaction in the period of the 
spectral flow. 

In this paper, we have been concerned with the cases 
of rational twist angles in unit of 2ir. We think that 
the present conclusion can be applied even for irrational 
cases, although we could not yet find the exact solution 
for such cases. 
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Note added in proof 



We have learned that the spectral flow for the XX Z 
model with short-range interaction was studied in F. C. 
Alcaraz, M. N. Barber and M. T. Batchelar, Ann. Phys. 
182 (1988) 280 . 



FIG. 1. Exact spectral flow for the system with N = 6, 
M — 2 and Q — 2: = cj>f = cj>± for (a) and <j> = 4>t , <t>i =0 for 
(b). Lower 40 levels are shown. In order to change the value of 
cj> — p/q, integer p is varied consecutively with q — 100 being 
fixed. We have drawn the diagram up to <f> = 2 for simplicity. 
When we discuss the flow beyond (j> — 2, we can use the pe- 
riodicity of the full spectrum, for example, <f> — 0, 1, 2 should 
read (j> — 2, 3, 4. Note that the spectral flow exhibits singular 
cusps, which come from the special dispersion relation (26) 
for the present model. See also [6] for detail. 
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